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Abstract 

A technique for the full wave characterization of micro- 
strip open end discontinuities fabricated on uniaxial aniso- 
tropic substrates using potential theory is presented. The 
substrate to be analyzed is enclosed in a cut-off waveguide, 
with the anisotropic axis aligned perpendicular to the air- 
dielectric interface. A full description of the sources on the 
microstrip line is included with edge conditions built in. 
Extention to other discontinuities is discussed. 


Introduction 

While there is extensive data available on the micro- 
wave characterization of a variety of microstrip discontinuities 
using both quasi-static [1-3] and full wave techniques [4-6], the 
characterization has been restricted to isotropic substrates. To 
date, there is no published data regarding microstrip disconti- 
nuities patterned on anisotropic substrates. Some very useful 
microwave substrates however, like sapphire, are anisotropic 
and so any discontinuity structures fabricated on them can 
not be properly characterized by the techniques that have 
been developed for isotropic dielectric substrates. 

A technique for the full wave characterization of micro- 
strip open ends fabricated on lossless uniaxial anisotropic 
substrates has been developed and is presented here. It is 
based on a dynamic source reversal technique that uses poten- • 
tial theory [7], which is a generalization of the charge resversal 
technique introduced several years ago [1], The discontinuity 
is enclosed in a waveguide of infinite extent whose dimensions 
are such that the guide is cut-off for the propagating fre- 
quency on the microstrip. All sources on the microstrip are 
represented, and the technique does not require a model for 
the source exitation. 


Dynamic Source Reversal Technique 

The anisotropic axis of the substrate is aligned perpen- 
dicular to the air-dielectric interface as shown in Fig. 1. The 
anisotropic dielectric may be represented as a tensor quantity 
given by 

*(y) = *(y)i + [«y(y) - *(y)lw (0 

where I is the unit dyad and #e(y) = * y (y) = 1 for y > h. The 
microstrip line is assumed to be infinitely thin and located at 
a height y = h + . In terms of the sources on the microstrip 
line, the scalar and vector potentials, $ and A respectively, for 
the dielectric loaded waveguide may be determined from 

(V 2 + Kk,, 2 )^ = -I1 0 J x (2) 

(V 2 + <ck 0 2 )A, = - Mo J x (3) 

(V 2 + * y k 0 X = - Wo {k - l)*(h)«(y - h) 

+ j“M 0 « 0 (*y - K)9t/dy (4) 

[ic^/Ac 2 + d 2 /^ 2 ) + 9(*( y)B/dy)/dy + K 2 (y)k 0 2 ]$ 

= “ P/ ( o + M*, “ l)AyOWy - h) - jw(/c y - li^Aj/dy (5) 

The potentials appearing in Eqs. (2) to (5) are obtained 
from the appropriate Green’s Functions and the corresponding 
sources using 

A, J (x,h^) = /i 0 j j G X I (x,h,z;x',h,z')J X I (x',h, 2 ')dx'dz' (6) 
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e 0 *(x,h,z) = f J G^(x,h,z;x',h,z')/>(x',h,z')dx'dz' (7) 

From these potentials the electric field components are found 
using 

E x = — jo;A x - dfc/dx 
E y - -jo; Ay - db/dy 

E s = — jwAj - &/dz (8) 

where it is required that for this particular geometry E x and 
E z vanish on the microstrip. The fields thus obtained are 
expressed in terms of LSE and LSM modes of the dielectric 
loaded waveguide. The anisotropic effect appears only in the 
LSM mode terms, which are present in A and The LSE 
modes are unchanged from those obtained for the isotropic 
case. 

A complete set of dominant mode sources on the micro- 
strip are represented; the longitudinal and transverse currents, 
as well as the charge on the microstrip line, with appropriate 
edge conditions built in. For a wide range of practical open 
end discontinuities, a valid approximation is that J x = 0, so 
therefore A x = 0. As a result, only the boundary condition 
= 0 is required for this problem. 

A line terminated at z = 0, thus forming an open end, 
would create reflected dominant mode sources on the line, 
along with perturbed sources localized near the discontinuity. 
The total source distribution on an open end may be written 
as 

J,(x'A»0 = J o - Re +i/5 *') + Jx »(*'•*') (») 

p(x‘M = /> 0 ( x ')(e' j ^' + Re + j*') + Pl (x',z') (10) 

for z 9 < 0, where J Q x , p Q and ^ are the yet to be determined 
amplitudes and propagation constant, respectively of the 
dominant microstrip mode, R is an unknown reflection coeffi- 
cient, and Jj a and p^ represent the perturbed source ampli- 
tudes near the open end. Weighted Chebychev polynomials are 
used to represent the sources in x for both dominant and 
perturbed sources, while triangle and pulse functions are used 
to represent the perturbed sources in z. 

Equations (9) and (10) may be written as 

J^x'.h.z') = [j(l + R)J 0 I (x')(B in cos(/9z') - sin(/3z')) 

+ jJ 11 (x',z')] (l-U(z')l (11) 

p(x'4i,z') = 1(1 + RVotx'XBfennf/te') + cos(^z')) 

+ /»,(x',z')] • [1 - U(z')] (12) 

where U(z') is the Heavyside unit step function which is 0 
for z ' < 0 or 1 for z' > 0, and jB^ = (1 — R)/(l + R) is the 
normalized input suscept&nce for the open end. In Eqs. (ll) 


and (12) the dominant mode sources are assumed to exist for 
— oo < z' < oo. Then dominant mode sources for z' > 0 are 
subtracted away to create the terms that multiply the (1 + R) 
coefficient in (ll) and (12). Since the amplitudes of J t % and 
p l are arbitrary at this point, they may be defined so as to 
include the (l + R) term. Now the (1 + R) term is common 
to all of the source terms, so it may be normalized to 1.0. 
Using Eqs. (ll) and (12) in Eqs. (6) and (7) and substituting 
into Eq. (8) gives the electric field in terms of the sources on 
the open end microstrip line. When the requirement that E s = 
0 on the microstrip is enforced, the terms corresponding to the 
dominant mode on the infinite line already satisfy the bound- 
ary condition on the strip, so they drop out. The sources 
existing for z' > 0 may be considered “source reversed” terms 
which produce an impressed Geld in the reigon z < 0 but 
localized near the discontinuity. The factor (1 + R) can be 
included into the arbitrary amplitude of the dominant mode 
sources. Thus, apart from the unknown parameter B^, the 
dominant mode sources in z' > 0 produce a known forcing 
function in the strip for z' < 0. The electric Geld produced by 
the perturbed sources Jj and p v must cancel the tangential 
component of the applied Geld for z < 0. A modiGed perturba- 
tion technique [8] is used to determine the unknown dominant 
mode amplitudes and propagation constant for an inGnite line. 
The method of moments is then used to reduce the resulting 
integral equation to a matrix equation which can be solved for 
the unknown input admittance B^, as well as for J x % and p v 
Only one matrix inversion is required to Gnd B in . 

Results 

The pulse width, A, of the expansion functions was 
chosen to be 0.32 mm at f = 2.0 GHz (or A ^ 0.0053 X for 
sapphire). This pulse width guaranteed a converged value for 
®in P] f° r examples presented. To verify the accura- 

cy of the theory as well as the resulting program, the program 
was checked for the isotropic case, and was able to duplicate 
data obtained in [l] and [7]. 

Table I shows the results obtained using this technique 
for several different anisotropic substrates as a function of 
microstrip line width. The open circuit capacitance, is 
found using C K = B^/wZq, where the characteristic imped- 
ance Z 0 is obtained from a computer program developed on 
the basis of the theory presented in [8] for the characterization 
of inGnite microstrip lines with sidewalls, but no top cover. 

To justify using values of Z Q thus obtained, calculations per- 
formed for a microstrip line with an air dielectric showed less 
than 2% difference in Z Q values obtained with and without a 
top cover. Table II shows the variation of as a function of 
line width for sapphire, and compares the results obtained for 
a substrate with an isotropic dielectric constant of 9.4, as well 
as for an isotropic dielectric constant of 11.6. Table III shows 
the effects of Gxed waveguide dimensions on B in and as a 
function of frequency of the propagating microstrip mode. For 
low frequencies, B^ varies linearly with frequency, starting to 
deviate as the frequency increases, this effect becomes more 
pronounced until the cut off frequency of the E^j waveguide 
mode is reached. This effect may be overcome by either fre- 
quency scaling the input parameters or by adjusting the 
waveguide dimensions accordingly. 
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The BASIC computer program developed to implement 
this technique can be executed on a personal computer with as 
little as 640K of RAM. Other discontinuity structures can be 
characterized in a similar manner. The technique is computa- 
tionally efficient, there is no need to model the source excita- 
tion, and the admittance can be solved for directly in the case 
of a one port network. All integrals involving the expansion 
and testing functions are performed analytically so no numeri- 
cal integrations are necessary, and the dominant portions of 
slowly converging series can be extracted and summed into 
closed form. 

This technique can be extended to rapidly and accu- 
rately characterize a number of other commonly used disconti- 
nuity structures, especially “coaxial” two port structures such 
as asymmetrical gaps and steps in width. To characterize a 
two port structure in terms of an equivalent “Tee” or “Pi” 
network, the Tangent Plane method [9] can be used to extract 
parameter values. 
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TABLE I.— OPEN CIRCUIT CAPACITANCE, C^, FOR SEVERAL ANISO- 
TROPIC DIELECTRIC MATERIALS. FREQUENCY = 2.0 GHz, h = 1.0 mm, 
b = 11 mm, 2a = 20 mm FOR W/h < 4.0 ELSE 2a = 100V/h).Z() VALUES 
OBTAINED FROM REF. 8. UNITS ARE pF/METER FOR C^. 



TABLE 0.— VARIATION OF C^ AS A FUNCTION OF LINE 
WIDTH FOR SAPPHIRE, k = 9.4, k y = 1 1 .6 COMPARED TO 
THAT OF AN ISOTROPIC DIELECTRIC. FREQUENCY 
= 2.0 GHz, h = 1.0 mm, b = 1 1 mm, 2a = 20 mm FOR 
W/h < 4.0, ELSE 2a = 10(W/h). UNITS ARE pF/METER FOR C^. 


W/h 


0J25 


03 


1.0 


2.0 


4.0 


6.0 


TABLE m.— VARIATION OF THE 
NORMALIZED INPUT SUSCEPT- 
ANCE, B in , AND OPEN CIRCUIT 
CAPACITANCE, C^, WITH 
FREQUENCY. WAVEGUIDE 
DIMENSIONS ARE 2a = 20 mm, 
b = 1 1 mm AND h = 1 mm. 


Frequency, 


C«(pF/m) 

GHz 
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